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One-sentence summary: The regions of the globe where soil moisture variations can 
affect precipitation are determined, for the first time, by a superior multi- model approach. 
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Abstract: Previous estimates of land-atmosphere interaction (the impact of soil moisture 
on precipitation) have been limited by a severe paucity of relevant observational data 
and by the model-dependence of the various computational estimates. To counter this 
limitation , a dozen climate modeling groups have recently performed the same highly- 
controlled numerical experiment as part of a coordinated intercomparison project. This 
allows, for the first time ever, a superior multi-model approach to the estimation of the 
regions on the globe where precipitation is affected by soil moisture anomalies during 
Northern Hemisphere summer. Such estimation has many potential benefits; it can 
contribute, for example, to seasonal rainfall prediction efforts. 


Atmospheric chaos severely limits the predictability of precipitation on seasonal 
timescales. Indeed, weather forecasts, which rely heavily on atmospheric initialization, 
rarely demonstrate skill beyond about a week. Hope for accurate seasonal forecasts lies 
with simulating atmospheric response to die slowly varying states of the ocean and land 
surface - components of the Earth system that can be predicted weeks to months in 
advance. A systematic response of the atmosphere to these boundary components \*ould 
naturally contribute skill to seasonal prediction. 

The critical importance of the ocean surface in this regard is well known 1 . Ocean 
temperature anomalies can be predicted a year or more in advance 2 . Furthermore, the 
atmosphere responds particularly strongly (and predictably) to ocean temperature 
anomalies in certain regions - “hot spots” of ocean-atmosphere coupling. The eastern 


equatorial Pacific is the most famous oceanic hot spot, playing a key role in the El Nino - 
La Nina cycle 3 . 

Another potentially useful slow component of the Earth system is soil moisture, which 
can influence weather through its impact on evaporation and other surface energy fluxes. 
Soil moisture anomalies can persist for months 4 , and though a paucity of observations 
prevents an unambiguous demonstration of soil moisture impacts on precipitation 5 , such 
impacts are often seen in atmospheric general circulation model (AGCM) studies 6 1 . In 
fact, some AGCM studies suggest that in continental midlatitudes during summer, 
oceanic impacts on precipitation are small relative to soil moisture impacts 8 . 

This suggests a question: Are there specific locations on the Earth’s surface for which 
soil moisture anomalies have an especially significant impact on precipitation? The 
identification of such hot spots would have tremendous implications for the design of 
seasonal prediction systems and for the associated development of ground-based and 
satellite-based soil moisture monitoring strategies, if such impacts are found to be local 
In a broader sense, such identification is clearly critical for understanding Earth’s climate 
system and the limits of predictability therein. 

While AGCM studies 9 10 11 12 and even numerical weather prediction model studies 13 
have addressed this question, published results are based on different experimental 
designs and reflect distinctive features of different model parameterizations. The 
coupling question, however, was recently addressed en masse by a dozen AGCM 



groups 14 , all performing the same highly-controlled numerical experiment. The 
experiments were coordinated by GLACE, the Global Land- Atmosphere Coupling 
Experiment 15 . Each model contributing to GLACE generated several ensembles of 
boreal summer (June through August) simulations designed to quantify that model’s 
land-atmosphere coupling strength 16 for that season. By combining the results across 
these models, we eliminate much of the undesired individual model dependence. We 
obtain, in effect, a unique result: a multi- model average depiction of the global 
distribution of land-atmosphere coupling strength. Given the limitations of the 
observational data, both now and in the foreseeable future, such a multi- model estimate 
of coupling strength distribution is arguably the best estimate attainable. 

Each GLACE participant performed an ensemble of sixteen simulations in which soil 
moisture varied between the simulations, and another ensemble in which the 
geographically- varying time series of subsurface soil moisture was forced to be the same 
across the sixteen simulations 1 1 . Coupling strength- the degree to which all prescribed 
boundary conditions affect some atmospheric quantity X- can be estimated 18 for each of 
the two ensembles with the diagnostic Q : 

f2 = (16CT 2 <x>-a 2 x ) / 15 c^x , 

where c^x is the intra-ensemble variance of X and o 2 <x> is the corresponding variance of 
the ensemble- mean time series - the single time series generated by averaging across the 
16 ensemble members at each time interval, chosen here to be six days. We are 



Figure 1. Mean of the land-atmosphere coupling strength diagnostic (the Q difference, 
dimensionless) for boreal summer across models participating in GLACE. 


The positions of the hot spots are not unexpected 8 19 , particularly if the soil moisture 
influence is presumed to be local rather than remote. In wet climates, for which soil 
water is plentiful, evaporation is controlled not by soil moisture but by net radiative 
energy. Thus, in wet climates, precipitation should be insensitive to local soil moisture 
variations. In dry climates, on the other hand, evaporation rates are sensitive to soil 
moisture but are generally too small to affect precipitation generation; in any case, the 
atmosphere in dry regions is predisposed to limit precipitation Only in the transition 
zones between wet and dry climates, where the atmosphere is amenable to precipitation 
generation (in particular, where boundary layer moisture can trigger moist convection? 0 ) 
and where evaporation is suitably high but still sensitive to soil moisture, can we expect 
soil moisture to influence precipitation The major hot spots shown lie mainly in such 


transition zones. 


interested, of course, in precipitation; to reduce noise, though, we take X to be the natural 
logarithm of the precipitation. Performing statistics on the logarithms of precipitation is a 
common practice in hydrology and meteorology, since unmodified precipitation 
distributions tend to be highly skewed. 

A study of the equation shows that outside of sampling error, £2 should vary from 0 to 1, 
with higher values implying a higher impact of the atmosphere’s boundary conditions on 
precipitation. To isolate soil moisture’s impact on precipitation from that of the other 
forcings, such as time-varying ocean temperatures, we compute the difference in the £2 
values between the two ensembles. In simple terms, this £2 difference describes the 
fraction of the precipitation variance explained by variations in soil moisture. 

Figure 1 shows the global map of the coupling strength diagnostic - the £2 difference - 
averaged across all of the participating models in GLACE Several hot spots appear, 
most notably the central Great Plains of North America, the Sahel, equatorial Africa, and 
India. Less intense hot spots appear in South America, central Asia, and China. 
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One must not assume from the map that all of the GLACE models place hot spots in the 
regions indicated. There exists extensive inter- model variability in the strength and 
positioning of the hot spots, a reflection of ongoing uncertainty in the proper way to 
represent the physical processes defining land-atmosphere coupling strength. Still, 
despite this variability, certain patterns appear frequently. These are the patterns that 
survive the averaging process and appear in Figure 1 . Each hot spot shown in the figure 
was produced, with varying magnitudes, by several of the participating models. 

The hot spots indicate where a global initialization of soil moisture may enhance 
precipitation prediction skill, at least during northern hemisphere summer 21 . Under the 
assumption that the soil moisture impacts are predominantly local, die hot spots indicate 
where tire routine monitoring of soil moisture, with both ground-based and space-based 
systems, will yield the greatest return in boreal summer seasonal forecasting. The hot 
spots are, in a sense, land-surface analogues to the ocean’s “El Nino hot spot” in the 
eastern tropical Pacific. 
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Popular Summary 

"Hotspots of Land Atmosphere Coupling”, by R. Koster, P. Dirmeyer, Z. Guo, and 21 
others. (Submitted to Science ) 


A wetter-than-usual soil may lead to higher-than usual evaporation, which in turn may 
lead to increased precipitation. This soil moisture - precipitation connection, if verified 
and properly utilized, would contribute significantly to seasonal forecasting efforts. 
Seasonal forecasters could then take advantage of the fact that initialized soil moisture 
anomalies can persist for months. 

The problem with verifying the soil moisture - precipitation connection with 
observational data is that the required data on the large scale do not exist and are 
logistically impossible to obtain. Climatologists have thus relied instead on modeling 
studies to quantify the connection. These modeling studies have their own limitations, 
however, most notably, the results can be strongly model dependent. 

GLACE is an international intercomparison project designed to quantify the strength of 
the soil moisture - precipitation connection (the "coupling strength") across a broad 
range of atmospheric general circulation models. Through GLACE, we find that the 
different models do indeed show a broad disparity in coupling strength distribution. 
GLACE, however, also provides an intriguing result. Despite the intermodel disparity, 
certain areas of the Earth show a large coupling strength in many models, suggesting that 
the existence of significant coupling strength in these areas is not so model-dependent. 
Given the lack of observational data, such a multi-model determination of areas with 
strong coupling strength is arguably the best estimate of such areas attainable by any 
method. 

To highlight these areas objectively, we average the coupling strength distributions across 
the twelve models participating in GLACE. The areas that survive the averaging process 
— i.e., the areas that are arguably not model-dependent — lie in the Great Plains of the 
United States, in the Sahel, and in northern India. Additional regions with lower 
coupling strength appear in central Asia and just south of the Amazon. If the soil 
moisture - precipitation connection in nature is mostly local, these "hotspots" of coupling 
show where soil moisture measurement with in situ or satellite sensors would provide the 
most benefit to seasonal forecasts. 


